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CN ' Abstract 

Recent progress in Heavy Flavour Production phenomenology are dis- 
cussed. In particular, the long-standing discrepancy of the Tevatron b 
production data is considered. It is shown that a better use of e + e~ 
data in constraining the effect of fragmentation reduces the discrepancy 
considerably. 
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The phenomenology of heavy flavour production has been so far a mixed success. 



(nJ ! 1 Introduction 

CO 

I Although the order of magnitude of the total cross sections, and the shape of 

CNI ' differential distribution is reasonably predicted, large discrepancy are present, 

especially for b production. 

The theoretical framework for the description of heavy flavour production is 
the QCD improved parton model. Besides the well-established NLO corrections 
to the inclusive production of heavy quark in hadron-hadron |l], |^, ^|, ^] , hadron- 
r"| , photon g, Q], and photon-photon collisions much theoretical work has 

been done in the resummation of contributions enhanced in certain regions of 
phase space: the Sudakov region, the large transverse momentum region and 



X 

5^ ■ the small- a; region. 



The theoretical effort involved is justified by the large variety of applications 
that heavy quark production physics has, in top, bottom and charm produc- 
tion. Besides the need of modeling these processes, heavy quark production is an 
important benchmark for testing QCD and parton model ideas, due to the rel- 
ative complexity of the production process, the large range of masses available, 
and the existence of different production environments, like e + e~ annihilation, 
hadron, photon-hadron and photon-photon collisions. 



* Invited talk at the 6-th International Workshop on Heavy Quarks and Leptons, May 
27/June 1, 2002, Vietri s/m, Salerno, Italy. 
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2 Total cross section for top and bottom 



Top production |j [10| has been a most remarkable success of the theoretical 
model. The measured cross section has been found in good agreement with 
theoretical calculations, as shown in fig. |]. Resummation of soft gluon effects 
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Figure 1: Results on top cross sections at the Tevatron. 

IPJ reduces the theoretical uncertainty in the cross section, pushing it toward 
the high side of the theoretical band. It remains, however, inside the theoretical 
band of the fixed order calculation, thus showing consistency with the estimated 
error. 

Recently, the HERA-B experiment has measured the bb total cross section 
[ p"2[ . This experiment is sensitive to the moderate transverse momentum region, 
where the bulk of the total cross section is concentrated. Since the production is 
(in a certain sense) close to threshold, resummation of Sudakov effects is useful 
also in this case. A comparison of the HERA-b measurement with a theoretical 
calculation is shown in fig. |[ The HERA-b result is compatible with the central 
value prediction, with the b pole mass around 4.75 GeV. Higher precision may 
constrain further the b quark mass. 
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Figure 2: bb total cross section versus theoretical predictions, including resum- 
mation of soft gluons. 

3 Differential distributions 

CDF has a longstanding disagreement with QCD in b production. A very recent 
publication of the B + differential cross section |l3[ has quantified the disagree- 
ment as a factor of 2.9 ± 0.2 ± 0.4 in the ratio of the measured cross section over 
the theoretical prediction. This discrepancy has been present since a long time, 
and it has been observed both in CDF and DO. Some authors JO] have argued 
that the discrepancy could be interpreted as a signal for Supersimmetry. 

Because of the large theoretical uncertainties, this discrepancy has been often 
downplayed. In fact, several effects may conspire to raise the b cross section 
to an appropriate value. It was early recognized that small- x effects may be 
important in b production at the Tevatron jjj . Resummation of these effects 
[ fH)[ p~G| leads to an increase of the cross section of the order of 30%. Threshold 
effects are small in this case PH[ , being below 15%, but they are nevertheless 
positive. Resummation of large logp T yields an increase in the p T spectrum 
of the order of 20% in the intermediate p T region jl7| . The presence of many 
possible enhancements should not lead, however, to excessive optimism. First 
of all, it is not clear whether these effect can be added up without overcounting. 
Furthermore, they arc all higher order effects, and thus should not push the 
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cross section too far out of the theoretical band, which includes estimates of 
unknown higher order effects. 

It has been observed since some time that an improper understanding of 
fragmentation effects may be one of the causes of the Tevatron discrepancy. 
This possibility stems from the fact that b quark cross sections are in reasonable 
agreement with the Tevatron measurements of the B meson spectrum, while 
the cross section obtained by applying a standard fragmentation function of the 
Peterson form with e = 0.006 to the quark cross section yields too soft a 
spectrum, as can be seen in fig. |^. In rcf. []l9| it was suggested to study b quark 
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Figure 3: Quark versus Hadron differential cross section compared to CDF data. 
The Hadron cross section is obtained by convoluting the quark cross section with 
a Peterson fragmentation function, with e = 0.006. 

jets rather than B meson's distributions. In fact, the jet momentum should 
be less sensitive to fragmentation effects than the hadron momentum. A DO 
study [^0| has demonstrated that by considering b jets instead of B hadrons the 
agreement between theory and data improves considerably (sec fig. ^). 

4 Fragmentation 

If fragmentation effects are so important, they should be assessed very carefully 
before claiming a discrepancy with data. Non-perturbative fragmentation is 
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Figure 4: DO b-jet data compared with the calculation of ref. U 



usually introduced by writing the hadron-lcvel cross section for B mesons as 

^ dp T dz-—D(z) S(p T — zp T ) , (1) 



dp T J dpi 

where p T is the hadron, and p T is the quark transverse momentum. Tradition- 
ally, D(z) is given by the Peterson form []18| 

z(l - z) 2 

D(z) oc , n [ \ 2 ' , 2 • (2) 
((1 - z) 2 + ez) 2 

Strictly speaking, this procedure is justifiable only in the limit of large transverse 
momenta (i.e. pr 3> nm,). It can however be consider a rough model of the 
effect of fragmentation at moderate and small transverse momenta, as long as 
the variable z is referred to the momentum or the kinetic energy, rather than 
the energy or the + component, of the quark momentum. 

Eq. (||) is often used with e = 0.006 for bottom and e = 0.06 for charm, in 
conjunction with shower Montccarlo programs. These values of e were obtained 
from fits to inclusive heavy flavour production distributions in e + e~ collisions. 
The shower Montecarlo accounts for hard fragmentation effects (i.e. hard gluon 
radiation), where hard means from the typical transverse momentum of the 
process down to some cutoff, of the order of the hadron mass. This accounts for 
scaling violation in the inclusive cross section for the production heavy quarksQ 

The inclusive cross section for single particle production in e + e _ collisions is sometimes 



5 



It has become common practice to adopt the value e = 0.006 also in con- 
junction with NLO calculations of b production cross sections. This practice is 
not totally justified, since the procedure used to obtain this value of e does not 
match the accuracy of the NLO calculation. A better approach is possible; one 
can use NLL resummation of transverse momentum logarithms in the context 
of e + e~ data in order to extract an appropriate non-perturbative fragmentation 
function, and supplement NLO heavy flavour production calculations with NLL 
resummation, in order to correctly account for scaling violation. The theoretical 
tools to do this arc in fact available in the literature: 

1 Single inclusive particle production in hadronic collisions pl[ . Single 
hadron production are described in term of NLO single parton cross sec- 
tion convoluted with a NLL fragmentation function; 

2 Heavy quark Fragmentation Function [p2f ; a method for the computation 
of the heavy quark fragmentation function at all orders in perturbation 
theory is developed, and applied at NLL. Several applications in e + e~ 
physics have appeared [^3], ^4], ^5|. p6fl . 

3 Single inclusive heavy quark production at large px [ p7{ ; item 1 and 2 are 
combined to give a NLL resummation of transverse momentum logarithms 
in heavy quark production; 

4 FONLL calculation of single inclusive heavy quark production; item 3 is 
merged without overcounting with standard NLO calculations. This pro- 
cedure has been implemented both in hadroproduction Jl7| and in photo- 
production |2^, |2^| . 

When using these ingredients, it was found that the CDF discrepancy is con- 
siderably reduced |30). It can be quantified as a factor of 1.7 ± 0.5 (theory) ± 
0.5 (expt). Reduction in the discrepancy is due to four basic points: 

1 FONLL calculation brings about a 20% increase in the intermediate px 
region. 

2 The fragmentation step for going from a "perturbative" b quark to a B 
hadron seems to be too strong at small px in the CDF implementation. 
We get a 20% increase in the small px region. 

3 A Peterson fragmentation function with e = 0.006 is too soft (present 
e + e~ data favour values around 0.002); this is a 20% effect. 

also called fragmentation function, which is often a source of confusion 
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Figure 5: Comparison of CDF data with the calculation of ref. Jsd[]. 

4 A more accurate use of e + e~ input data on the fragmentation function, 
described in the following section, brings about another 20% effect. 

The net effect is 1.2 3 = 1.7 and 2.9/1.7 = 1.7, so that the 2.9 CDF discrepancy 
becomes our 1.7. 

5 A better use of e + e data on b fragmentation 

D(z) is extracted from fits to e + e~ data. One has 



dp 



dp 



measured computed 
From the above equation it is easy to show that 



fitted 



= (xr 1 ) d 



N 



where 



x B 



Xb 



P 



-Djv = / dzz 



JV-l 



D(z). 



(4) 
(5) 



Pmax Pmax 

It turns out that for the computation of the hadronic cross sections, only the 
first few moments of the fragmentation function are important. This follows 
from the fact that the heavy quark hadroproduction cross section is a steeply 
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Figure 6: As in fig. ratios of data to theory. 
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falling function of the transverse momentum, following roughly a power law. 
Assuming the form 

da{H x H 2 -^b + X) 



dp T 



Ap T ' 



one finds 
da{HxH 2 -> b + X) 



dzdp T D(z) — S(p T — zp T ) 



A 



-D r , 



da 

dp T 



(6) 



(7) 



This simple fact was noticed long ago 31 , 32j ■ The exponent n ranges from 3 to 
5 in hadronic collisions. In ref. it was shown that eq. (0) is also an excellent 
approximation to the exact convolution. 

When fitting the non-perturbative fragmentation function from e + e~ data, 
the low moments relevant for hadroproduction are not well fitted. The Peterson 
form with e = 0.002 gives a good fit to the shape of the inclusive cross section, 
but not to the moments. For the purpose of predicting hadronic cross section it 
is instead better to make sure that the first few moments are accurately fitted. 
In fig. 0, it is shown that a fit to the second moment alone with the simple form 
x(l — x) 13 , yields a good fit to all moments in the interesting region (i.e. between 
3 and 5). We also notice that the second moment fit (N = 2 fit from now on), 
the e = 0.002 and the e = 0.006 fits differ by 20% around the N = 4 moment. 

The reason why x fits do not yield good fits to moments is easy to explain. 




Figure 7: The measured (x^ -1 ), compared with the perturbative NLL calcu- 
lation supplemented with different non-pcrturbative fragmentation forms. The 
solid line is obtained using a one-parameter form fitted to the second moment. 



The large x region in the inclusive production of heavy flavoured hadrons in 
e + e" annihilation is a very difficult region to model. It is in fact affected by 
large perturbative (i.e. Sudakov) and non-perturbative effects. It is unlikely 
that the simple forms used commonly for the non-perturbative fragmentation 
function can accurately model these complex effects. Thus, the large x region 
is usually excluded from the fits. On the other hand, this region is the most 
important one for moments. 

From fig. [?], another important fact becomes apparent. When using NLL 
calculations for the perturbative part of the fragmentation function, the im- 
portance of the non-perturbative effects on moments is greatly reduced. This 
can be taken as an indication of the fact that the heavy quark fragmentation 
function is, to a considerable extent, perturbatively calculable. 



6 Other B production data 



The procedure of ref. [|30| ma y affect other b production data at hadron colliders. 
Here I consider, as an example, the inclusive muon cross section in the central 
and forward region, where the muons come from the scmilcptonic decays of the b 
flavoured hadrons. These distributions were measured by the DO^] experiment, 
and found to be higher than theoretical prediction by up to a factor of 4 in the 
forward region. The result of the calculation performed with the method of 
ref. (30) is reported in fig. ||. As one can see, the agreement with data is now 
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Figure 8: DO forward muon data compared with a calculation performed using 
the method of ref. |?^/. 

quite remarkable. 

7 Conclusions 

The theory of heavy flavour production seems to give a good qualitative de- 
scription of the available data. In the case of top production, the comparison 
between theory and experiment is satisfactory also at a quantitative level. 

Recent progress has taken place in the field of b hadroproduction. The 
HERA-b experiment has provided a cross section for b production at relatively 
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low CM energy. Some progress in understanding the role of fragmentation has 
considerably reduced the longstanding problem of the b momentum spectrum 
at the Tevatron. 

Major problems do remain in the (perhaps less developed) areas of bottom 
production in 77 and 7p collisions. A discussion of these problems is given in 
ref. @. 
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